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A minimal model for multipole orders in CeBg shows that degeneracy of the quadrupole order param- 
eters and strong spin-orbit coupling lead to peculiar temperature and magnetic-field dependences of the 
X-ray reflection intensity at superlattice Bragg points. Furthermore, the intensity depends sensitively on 
the surface direction. These theoretical results explain naturally recent X-ray experiments in phases II and 
III of CeBg. It is predicted that under weak magnetic field perpendicular to the (111) surface, the reflection 
intensity should change non-monotonically as a function of temperature. 
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In CeBe and related systems, presence of both orbital 
and magnetic degrees of freedom brings about rich struc- 
tures in the phase diagram. In zero field, CeBe turns into 
the antiferro quadrupole (AFQ) ordered phase (phase II) at 
Jq = 3.4K and turns into antiferromagnetic (AFM) ordered 
phase (phase III) at - 2.3K.' The order parameter in 
phase II is the Fsg-type quadrupole moment. In phase III, the 
non-collinear magnetic structure is described by four wave 
numbers: ki = [1/4, 1/4, 1/2], = [1/4, -1/4, 1/2], Jtz = 
[1/4, 1/4,0], jtj = [1/4,-1/4,0]. 

Recently, X-ray scattering has been utilized as a powerful 
probe to detect orbital orderings.^ For CeBg, Yakhou et al. 
have performed resonant and non-resonant scattering experi- 
ments and found superlattice reflections in phases II and III.^ 
Nakao et al. have identified the boundary between phases I 
and II in a magnetic field by resonant scattering.'* More re- 
cently, Tanaka et al. have reported unexpected temperature 
and magnetic field dependences using non-resonant X-ray 
scattering.^' ^ Although the superlattice spots emerge below 
Tq, the intensity of (n/2, n/2, n/2) reflections with n odd in- 
tegers remains small in Tn < T < Tq. The intensity increases 
almost stepwise below T;^.^-^ Furthermore the small intensity 
in phase II is suppressed by application of a magnetic field of 
as small as 0. 1 T, while the suppression in phase III requires 
an order of magnitude larger magnetic field. ^ These features 
seem strange at first sight since the staggered quadrupole mo- 
ments probed by X-rays should already be present in phase II, 
and should not change significantly below the Neel tempera- 
ture Tn. In this paper we demonstrate how X-rays probe the 
coupling among dipole, quadrupole and octupole moments, 
and show that consideration of a quasi-continuous symmetry 
of the quadrupole order parameters provides a natural expla- 
nation of the experimental observations. 

In CeBg, a localized 4f electron of a Ce^^ ion has the quartet 
crystalline electric field (CEF) ground state, which is called Fg 
and is well below the excited CEF level Fy.^ An orbital pair 
in the Fg level are given by 



l+T> 



(1) 



in terms of eigenstates of J.. The Kramers partners of |± f) 
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are written as |± X), and are obtained by reversing the sign of 
J- in eq. (1). To describe multipole operators, we introduce 
two kinds of pseudo spin operators cr and r as 



r|±T> = ±l±T), T-^|±i) = ±|±i>, 

(r-^|±T> = l±T>, cr^|±i> = -|±i). 
Then the magnetic moment M is given by 



M = //B 2 (o"'- + ' 



(2) 



(3) 



where j/, describes the orbital dependent part, and is given by 
T] - {ri^cr\ rj'cr- , t'ct^), (4) 

with rf - (+ V3t^ - t")/2. The Fj^-type quadrupole moment 
has components 



T-'cr. 



(5) 



It is convenient to introduce a vector fi - (0,^, O^.x, Ow). The 
octupole moment with the F5,, symmetry has also three com- 
ponents given by 



(6) 



with f± = -(r' + V3t^)/2. 

We work with a RKKY type multipole Hamiltonian^ '" 
which reproduces phases II and III with minimum number of 
interactions. The model under magnetic field H is given by 



<ij) 

iU) 77' 



(7) 



where {ij) denotes a nearest neighbor pair, and [ij] denotes 
a next-nearest neighbor pair We have introduced the next- 
nearest interaction of the pseudo-dipole type: 



1^7 7 



(8) 



KAuiid''-'' -3<.n[.)/12, 

where n,j is the unit vector across the next-nearest neighbor- 
ing sites / and j. This interaction stabilizes phase III with the 
peculiar pattern of dipole moments. For simplicity, we do not 
consider such part of dipole interactions that comes from cr, 
in eq. (3), nor the F5„ type pseudo-dipole interaction, which 
is known to be important for realizing phase III' in magnetic 
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Table I. Typical directions of surface and magnetic field. 



configuration 


S(lll) 


S(llO) 


surface direction 


(1,1,1) 


(1,1,0) 


magnetic field direction 


(2,1,1) 


(1,1,0) 




(n/2,n/2,n/2) 


(/i/2,/i/2,l/2) 




n = 1,3,5,7,9,11 


ft = 5,7,9, 11 



field larger than 1 T.'" Hence discussion of phase III' is out 
of the scope of this paper. We assume that Dsg is positive and 
is the largest among all interactions. Then the Fs^-type AFQ 
order first sets in from the paramagnetic phase. The term 
is introduced to simulate the symmetry breaking due to the 
surface, as explained later. 

We study this Hamiltonian by the mean field theory with 
superlattice structures up to Vs X Vs X 2 supercell, which 
can describe the magnetic structure in phase III. We take the 
energy unit as the quadrupole ordering temperature: Tq = 1 
or, equivalently, Dsg - 1/6. With the choice D4„2 - Dsu = 
0.9D5g and K4u2 = 0.5, we obtain the magnetic ordering tem- 
perature Tn = 0.47 with zero field. Our mean field theory 
indeed stabilizes the pattern of the dipole moments proposed 
in ref . 1 . 

In phase II, this model realizes the Ts^-type order parame- 
ter with three components (Oy^, O^, Oxy) at each site. In the 
mean field theory for eq. (7), there is no preferred direction 
for without magnetic field. We call this situation a quasi- 
continuous symmetry. Because of the spin-orbit coupling, an 
external magnetic field tends to aUgn n\\H. This coupUng is 
apparent in i/ in eq. (3). In other words, dipole and octupole 
moments combine to give t] under the cubic symmetry.'- In 
phase III, on the other hand, simultaneous presence of dipole 
and octupole orders restrict the direction of ^ even without 
magnetic field. In real CeBe, there should be various sym- 
metry breaking perturbations to fix the direction of n even in 
phase n. We consider in particular such effect of surface that 
is simulated by the following term: 



(9) 



where e is the unit vector normal to the surface, and £(> 0) as- 
sumes that the surface prefers the wave function extended par- 
allel to it. For example, the (0,0,1) surface prefers Oxy to other 
components Oy^, O^x- In phase II where AFQ is present, the 
(0,0,1) surface disfavors the AFQ component {Oxy)AFQ of the 
order parameter. This situation is analogous to the Neel state 
in the antiferromagnetic Heisenberg model, where a mag- 
netic field disfavors the magnetic moment parallel to the field. 
Since we do not go into the details of the surface region, we 
take the summation over ; in eq. (9) for the whole system. 

We now consider the consequences of symmetry breaking 
by the surface and by magnetic field on X-ray reflection. We 
take two typical conditions as shown in Table I, which cor- 
respond to recent experimental configurations.^'^ To calculate 
the Thomson scattering intensity, we need the structure factor 
F{k) defined by 



F{K) = 2 {a\ exp(i*c • r) \a) , 



(10) 



where R„ specifies a Ce site, and the coordinate r of a 4f elec- 
tron is measured from the Ce site. The momentum transfer is 



given by k. The state labeled a has the statistical weight pa. 
Then the scattering cross section is given by 



dcr\ 
dflj. 



mc 



(11) 



where e denotes polarization of the incident beam and e' that 
of the scattered beam." 

We expand \a) in terms of the total angular momentum ba- 
sis \JM) with J = 5/2 as \a) - Y^m cm \JM). Then we obtain 
in eq. (10) 

{a\ exp(iK • r) \a) 

= ^ Cj^cm' ^ ^ {lm\ exp(iK • r) \lm') 

MM' m,m' m. 



X {JM\lmsms) [Im' sms\JM') , 



(12) 



where (JM\lmsm^) is a Clebsch-Gordan coefficient with / = 3 
and s = 1/2. The wave function (r|/m) is factorized into the 
radial and angular parts as {r\lm) - f(r)Y'^(f), where r = \r\ 
and f is the sohd angle. From this we obtain 

{lm\ exp(i*c • r) \lm') 



K 

X 2 i-^fY^aiK) 

Q=-K 

I K I \l I K I 
)1 -OT Q m! 



x(-l)'' 



(13) 



(14) 



Here we have introduced the notation (jK(i<)) by 

(Mk)) = J ArjK{Kr)f{ry, 

where jK{Kr) is a spherical Bessel function of order K. We use 
the data of {jKiK}) calculated by Freeman and Desclaux with 
Dirac-Fock method.'^ 

Before considering the effects of magnetic field and sur- 
face, we impose by hand an AFQ order with only Oxx and 
derive the X-ray scattering intensity. This artifice helps us to 
identify the relationship between the AFQ order and k. The 
wave functions diagonaUzing Oxy are written as 

1 



|AT> = -p(l+ T)+i|- T», 

V2 

IB T) = ^(1+ T>-i|- T)), 

V2 



(15) 



and their Kramers partners \A I) and \B I). The eigenvalue of 
Oxy is ±1 for \A ti), and +1 for \B ti>. We take the staggered 
(G-type) AFQ order with complete occupation of \A t) for the 
sublattice A, and that of \B t) for the sublattice B. Since the 
X-ray intensity does not depend on the spin direction, this fer- 
romagnetic configuration serves to analyze basic features of 
the intensity. Figure 1 shows the intensities of superlattice re- 
flections against sin0//l = K/{4n), where A is the wavelength 
of X-ray, and 26 is the scattering angle. The intensity is de- 
fined as the scattering cross section with |e' • e| = 1, and 
is the classical electron radius = e^/mc^. From Fig. 1, it 
can be seen that the intensities have strong dependence on the 
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Fig. 1. The superlattice reflection intensities in the O^y AFQ state with n 
odd integer against sin 9/ A. 



direction of The reason for this behavior is explained quali- 
tatively as follows: In the AFQ state, the structure factor with 
K being a superlattice vector is written as F{k) oc fj^iK)- fB{K), 
where fA(B)(K) denotes the form factor of the A{B) sites. There 
are high-symmetry directions from which A and B sublattices 
are seen as having the same projected charge density. If a scat- 
tering vector is along such direction, we obtain F(k) = since 
fA(K) = In the case of Oxy AFQ state, such directions 

are x and y axes. Therefore the (1/2, n/2, 1/2) intensity with 
large n is small in Fig. 1. On the other hand, we also have 
fAiK) = fBiic) for K II (0,0, 1) since the different distribution 
of charge density in the xy plane does not survive integration 
over this plane. Thus the (1 /2, 1 /2, n/2) intensity is also small 
with large n. 

To the contrary, the difference of the charge density should 
be seen most effectively from the direction parallel to (1, 1, 0). 
In our calculation shown in Fig. 1, (n/2, n/2, 1/2) reflection 
has a deep minimum around K/4n ~ 0.6A"'. This comes from 
interference of K - 2 and K = 4 contributions as discussed in 
ref.l3. Namely, we obtain 

F(k) oc k,Ky{j2iK) + 5l2{lkl - \)j4{K)], 

where ka is directional cosine along the a axis.'^' Because 
k^ becomes small as k increases in the (n/2, n/2, 1/2) direc- 
tion, the coefficient of 74 becomes negative and suppresses 
the A" = 2 contribution. Hence the scattering intensity be- 
comes minimum when these two contributions nearly cancel 
each other. 

We now proceed to mean-field treatment of phases II and 
III, and discuss X-ray intensities without assuming the AFQ 
pattern a priori. We tentatively take E - 5 X 10~^ in eq. (9). 
Figure 2 shows the temperature dependence of (n/2, n/2, n/2) 
reflections in the configuration S(lll). It can be seen that 
(n/2, n/2, n/2) reflections have no intensity in phase II. By 
the term simulating the surface effect, avoids the direction 
( 1 , 1 , 1 ) in phase II. This means that the order parameter does 
not have a component (9^ , + (9.^ + (9 „ , which alone contributes 
to a finite scattering intensity for (n/2, n/2, n/2) reflections. 
Other components lead to fAiic) - fsiK), and the intensity is 
zero even with a quadrupole order. 

On entering phase III, yu gains a component in the direc- 



Fig. 2. The temperature dependence of X-ray intensity for the configura- 
tion S(lll). The scattering vector is (n/2, n/2, n/2). The Intensity remains 
zero in phase II. 



tion of (1, 1, 1), since the surface term is no longer the only 
source to fix the direction of ^. In the domain where we have 
the dipole order of (1/4, +1/4, 1/2), the Oxy component of 
quadrupole moment is stabilized by the spin-orbit coupling. 
Another component Oy^. is stabilized in the (1/2, ±1/4, 1/4) 
domain, and O^x in the (1/4, 1/2, ±1/4) domain. These do- 
mains give same superlattice intensity in the configuration 
S(lll). In our mean-field theory, dipole and octupole order 
parameters develop continuously from zero below Tn. How- 
ever, the direction of fi changes discontinuously at Tn, and 
the intensity changes discontinuously. Actual experimental 
results show small but finite intensities of (n/2, n/2, n/2) in 
phase 11 with the configuration S(lll).^ We interpret this 
feature in terms of slight mixture of unfavorable component 
Oyz + Ozx + Oxy in the sample by imperfections other than 
surface and by multidomain effects. 

In the configuration S(l 10), the surface favors ± (1, 1, 0). 
The X-ray scattering with k \\ (1,1,0) probes the quadrupole 
component Oxy = Mz- Hence (i \\ (0, 0, 1), which is favored by 
eq. (9), contributes to the scattering. On the other hand, an- 
other favored component || (T, 1, 0) is not probed by X-rays 
with a: II (1, 1,0). In this paper, we take the simplest approach 
to take both quadrupole configurations into account. Namely, 
we choose /j \\ (1,1, V2) which has equal weights of (1, 1, 0) 
and (0,0, 1) configurations. Figure 3 shows the temperature 
dependence of (n/2, n/2, 1 /2) reflections for the configuration 
S(llO). The intensities of (n/2, n/2, 1/2) reflections become 
finite below Tq and increase as the temperature is decreased. 
On entering phase III, the intensity increases discontinuously 
by the same reason in the case of S(l 1 1). 

A magnetic field can also influence the direction of fi by 
the spin-orbit interaction. Figure 4 shows the dependence of 
X-ray scattering intensities on magnetic field in the configu- 
ration S(llO). The intensity of (7/2,7/2, 1/2) reflection has 
a finite value without field, but is suppressed by a magnetic 
field. This is because the magnetic field along (1, 1,0) rotates 
ft so as to be parallel to H. Since fi in this direction is perpen- 
dicular to the scattering direction (7/2, 7/2, 1/2), the intensity 
vanishes although the quadrupole order is present. 

In the actual result for the configuration S(l 1 1), a smaU in- 
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Fig. 3. The temperature dependence of X-ray intensity in the configuration 
S(llO). The scattering vector is (n/2,n/2, 1/2). 



This is because the magnetic field induces the order parame- 
ter Oy^ + 0„ + Oxy. Since the dipole order in phase III will 
rotate /x from this preferred direction, the intensity should de- 
crease on entering phase III. We further predict that an inten- 
sity minimum of the superlattice reflection as a function of k 
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Fig. 5. The temperature dependence of X-ray intensity in the configura- 
tion S(lll) with H II (1, 1, 1) and hbH = O.ITq. The scattermg vector is 
(n/2,n/2,n/2). 



should be observed if (n/2, n/2, 1 /2) reflections are measured 
in phase III with the (1, 1, 0) surface. 
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Fig. 4. Temperature and magnetic field dependence of the X-ray intensity 
in the configuration S(llO). The scattering vector is (7/2,7/2, 1/2) and 
magnetic field is along (1,1, 0). 



tensity which is already present in phase II is suppressed by a 
small magnetic field. ^ This is explained in terms of rotation of 
fi toward H parallel to the surface. If the present mechanism 
is realistic, magnetic field along (I, I, I) should lead to reduc- 
tion of superlattice intensity on entering phase III. The result 
of model calculation is shown in Fig. 5. 

To summarize, we have clarified the efl'ects of quasi- 
continuous symmetry of the quadrupole moments on X-ray 
scattering from CeBs. The coupUngs with dipole and octupole 
moments lead to unexpected features; large increase of the 
superlattice intensity on entering phase III, and strong sup- 
pression of intensity by magnetic field. Our theory provides 
an interpretation of recent experimental results.^- ^ If magnetic 
field is applied perpendicular to the (1,1,1) surface, our the- 
ory predicts increase of the scattering intensity in phase II. 
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